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We give a general procedure to obtain non perturbative evolution operators in closed form for 
quantized linearly polarized two Killing vector reductions of general relativity with a cosmological 
interpretation. We study the representation of these operators in Fock spaces and discuss in detail 
the conditions leading to unitary evolutions. 
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^ i I. INTRODUCTION 

The study of symmetry reductions has been a fruitful way to gain valuable insights into the behavior of general 
relativity in its quantum regime. This is so because they provide tractable models where it is possible to make 
computations and obtain concrete predictions, at least in a qualitative way, about relevant features of a full blown 
^ ■ theory of quantum gravity. Two Killing vector reductions are specially appealing because they have local degrees of 
freedom and (restricted) diffeomorphism invariance, two of the features of general relativity that lie at the heart of the 
difficulties encountered in its quantization. These reductions differ from each other in the spacetime topology and the 
action of the isometry group. When the Killing fields commute and are hypersurface orthogonal the models become 
specially simple and it is possible to solve them in closed form in a straightforward way. In some cases the reduced 
model describes the propagation of linearly polarized wave-like modes in a spacetime with non-compact spatial slices, 
the so called Einstein-Rosen waves [1]. Here the symmetry group is 1 x U(l) and the spacetime is topologically K 4 . 
In other situations, when the symmetry group is U(l) x 1/(1), it provides cosmological models with initial singularities 
and an assortment of spatial topologies: There are non-compact examples such as the model introduced in [2] (referred 
to in the following as the Schmidt model) and compact ones such as the well known Gowdy T 3 model [3, 4], among 
many others. All these symmetric sectors of vacuum general relativity share the fact that the reduced phase space 
• • . can be parametrized by a scalar field in 1+1 dimensions and its canonically conjugate momentum (plus, eventually, 
some particle-like degrees of freedom). The main difference between them, as far as the quantization is concerned, is 
that the Hamiltonian for Einstein-Rosen waves [5] is time independent and has its origin in the surface terms needed 
to have a well-defined action principle whereas the Hamiltonian in the cosmological models is time dependent and is 
obtained through deparametrization [6, 7]. 

The quantization of linear Einstein-Rosen waves (free [5, 8, 9, 10] or coupled to matter [11]) is fairly well understood. 
The quantum unitary evolution operator can be obtained in closed form in a straightforward way and can be used for 
a number of purposes leading to physical applications such as the discussion of the existence of large quantum gravity 
effects [10] or the study of the microcausality of the system [12]. Models of the cosmological type, specially the Gowdy 
T 3 one, have been much harder to crack even though they received a lot of attention already in the seventies by Misner 
[13], Berger [14], and other authors. Although the first attempts to quantize the system were largely successful [15] 
and key technical insights were introduced already at these stage some features of the formalism were not completely 
satisfactory. In particular, the reliance on a Hilbert space built as a tensor product of infinite, one-particle Hilbert 
spaces is problematic since it is known that such a tensor product is not separable and the representation of the 
canonical commutation relations that it provides is reducible [16]. A renewed interest in the quantization of the 
Einstein-Rosen, Gowdy, and Schmidt models arose in the nineties [5, 6, 7, 9, 10] when, instead of the Dirac approach 
previously used, the quantization of the system was done by gauge fixing. At the time it was thought that they could 
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provide suitable testing grounds for loop quantum gravity and address some general issues related to the quantization 
of general relativity. However the finding that there were problems with the unitary implementation of dynamics in 
the Gowdy T 3 model (and also in the Schmidt case) [17, 18] was somewhat of a shock and was perceived as a potential 
drawback to its use as a toy model for quantum gravity. The situation has recently improved after Corichi et. al. 
[19, 20] have shown the existence of unitary evolution for a field closely related to the scalar that usually encodes the 
local gravitational degrees of freedom. 

The purpose of the present paper in this context is to obtain evolution operators explicitly in closed form written 
in terms of the basic objects, i.e. the field and momentum operators. This is done before choosing a specific Hilbcrt 
space representation for them. So, even if they are formal at this stage, they offer the possibility to explore different 
choices for the quantization of the system and, in particular, non-Fock representations (when available). Owing to 
this fact it is important to notice that our approach differs from the usual ones that make use of the Fock space -built 
from the solutions to the field equations- and a choice of complex structure to select the one-particle Hilbert space. 
As an application of our formalism we will discuss the unitary implemcntability of the time evolution of the system as 
done in [19] by changing the basic fields used to encode the physical degrees of freedom of the model. This provides 
an alternative point of view on this problem. 

The outline of the paper is the following. After this introduction we review in section II the Hamiltonian formalism 
for the U(l) x U(l) symmetric models considered in the paper; specifically we will discuss their derivation from an 
action principle, gauge fixing, and deparametrization. After this we will devote section III to obtain formal unitary 
evolution operators U(t,t ), defined in terms of abstract field and momentum. This will be done in a unified way 
for a family of time dependent Hamiltonians that include the Schmidt and Gowdy models. Under the condition that 
the above mentioned formal operators can be made unitary by choosing a suitable representation for the field on 
a Hilbert space they can be shown to satisfy the evolution equation ihd t U(t,t n ) = H(t)U(t,to) where H(t) is the 
time dependent Hamiltonian of the system. The technical details of the construction -that relies upon known results 
concerning the quantization of the time-dependent harmonic oscillator- will be left to appendix A. 

As an application of our scheme we discuss in section IV the problem of finding suitable Fock space representations 
for the field and its momentum that give rise to unitary evolution. At this stage the use of the auxiliary field considered 
in [14] and reintroduced in [19, 20] will play an important role that will be clarified within the present framework. 
We end the paper with our conclusions in section V and three appendices where we provide mathematical details on 
the obtention of the evolution operators and related topics. 

II. LINEARLY POLARIZED [7(1) x (7(1) MODELS: HAMILTONIAN FORMULATION 

The Gowdy and Schmidt models have been extensively used to study issues in quantum cosmology. In this context 
they are interesting because they share with the full theory of general relativity some of the features that renders its 
quantization highly non-trivial, in particular, they are genuine field theories with (some) diffeomorphism invariance. 
They have also been used to discuss physical issues such as the quantum behavior of the initial singularity in a 
non-homogeneous setting [14, 21]. 

A. The Midi-Superspace 

We will start by reviewing here some relevant facts about the Gowdy and Schmidt models, in particular, the 
symmetry reduction process and the deparametrization needed to define a convenient time-dependent Hamiltonian. 
The spacetimes (M, (4> g a b) of the Schmidt and Gowdy T 3 models [2, 4] are vacuum solutions to the Einstein equations 
characterized by their isometry group, in this case U(l) x 17(1). Once the group action is properly defined (smooth, 
effective, proper, and with no fixed points 1 ) on the connected and oriented spatial sections E 3 of M = K x S 3 the 
spatial topology is fixed to be of the form S 3 = X x S 1 x S 1 , where X = R for Schmidt and X = S 1 for Gowdy [22, 23]. 

We will focus here on the linearly polarized case where the isometry group is generated by a pair of mutually 
orthogonal, commuting, spacelike, and globally defined hypersurf ace- orthogonal Killing vector fields (£ a ,cr a ). The 
symmetry of the system implies that the Lie derivatives of the metric vanish: L^ (4) g a b = L a (l) g a b = 0. 

With these assumptions it is possible to exactly solve the vacuum Einstein equation 

W Rab = 0, (1) 



This condition rules out the S 1 X S 2 , S 3 , and lens space topologies for the Gowdy model. 
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where {i) R a b is the Ricci tensor associated with the Levi-Civita connection (4) V compatible with the 4-metric {i) g a b- 
Since the Killing fields are hypersurface orthogonal, the space of orbits M = R x X x S 1 defined by one of them, say 
£ a , can be identified as an embedded hypersurface in M. which is everywhere orthogonal to the (closed) orbits of £ a . 
The induced 3-metric {3) g a b on M can be written in terms of £ a := {4) g a b£, b and A := {4> g a b £ a £, b > as 

{3) g ab = (4> g ab - X-'Ub ■ (2) 

Let <3) i? a fc, <3) V Q , and (3) D be, respectively, the Ricci tensor, the metric connection, and the d'Alembert operator 
associated to (3> g a b- The equation (1) is then equivalent to the system [24] 

< 3 >DA = ^< 3 V b (< 3 >V a A) ( (3) V b A), 

(3) 

(3) Rab = 2i (3) V a (3) V fc A - ^ (< 3 >V a A) (< 3 »V b A) . 

The remaining Killing vector field a a satisfies L a (3) g a b = 0, L a \ = 0. 

It is well known that in order to further simplify the equations (3) we can use the conformal transformation 

gab ■= A <3) 5ab, (4) 

that is well defined since £ a never vanishes. Let R a b, V Q , and □ be, respectively, the Ricci tensor, the metric connection 
and the d'Alembert operator associated to the new 3-metric g a b- After the conformal rescaling (4), the system (3) is 
equivalent to 

□0 = 0, 

(5) 

Rah = (V„0) (V 6 0) , 

where we have defined the field V2(j) := log A. Therefore, the symmetry reduced models considered here can be 
thought of as (2+l)-general relativity coupled to a massless scalar field 0, with the additional symmetry defined by 
the Killing field a a . They can be derived from the 2+1 Einstein- Hilbert action for gravity coupled to a massless scalar 



that we will use as the starting point for the Hamiltonian formalism. The last term must be introduced in the non- 
compact case in order to have a well defined variational principle and involves an integration over the 2-dimensional 
asymptotic boundary of M. Here R is the scalar curvature of g ab , h and K are, respectively, the determinant of the 
induced metric and the trace of the second fundamental form of the boundary dM, and G3 is the Newton constant 
per unit length in the direction of the symmetry orbits. 

In order to respect the symmetry in the Hamiltonian analysis of (6) it is important to notice that, owing to the 
hypersurface orthogonality condition, the metric can be decomposed as 

gab = K b + T~ 2 (T a (T b , (7) 

where a a := g a b<J h , t 2 := <j a a a is the area density of the symmetry group orbits, and h ao denotes the induced metric 
of signature ( — h) on the 2-dimcnsional manifolds topologically R x X- that are everywhere orthogonal to the closed 
orbits of a a . The symmetry present in this case implies that the gradient of r is always timelike. We introduce now a 
foliation of R x X defined by spacelike level hypersurfaces of a suitable scalar function t. We also define a dynamical 
vector field t a such that t a \7 a t = 1. By using the unit, future-pointing, timclike, normal vector field to the foliation 
n a and the unit vector field x a compatible with the chosen orientation and tangent to each slice (topologically X) it 
is possible to write t a = Nn a + N x x a 1 where N and N x are the lapse and shift functions. The 2-metric h ab can be 
written then as h a b = -~n a rib + x a Xb, where n a and x a are respectively g a bn b and g a bX b - We introduce the additional 
non-unit vector field 

x a = e 7/2£a 

where 7 is an extra field, and impose the vanishing of the Lie brackets of the set of 

vector fields (t a ,x a ,a a ): 



d a N = d a N x = d al = 0, 
[a,x] a = [a,n] a = 0, 

Nei/ 2 [x, n} a = -n a d x N - x a (d x N x - dte 1 ' 2 ) 



(8) 
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By doing this it is possible to construct a global coordinate system for M: (t,x,a) e (0, oo) x X x S 1 . In these 
coordinates 

gab = (~N 2 + {N x ) 2 )\7 a t\7 b t + 2c~</ 2 N x \7 (a t\7 b) x + <?V a xV b x + t 2 V q ctV 6( t, (9) 

where the scalar fields N, N x , and 7 depend only on t and x. Thus, the midi-superspace under consideration consists 
of five real-valued smooth functions (N,N X ,T,^/,(f) that depend only on (t,x) and satisfy the Einstein-Klein-Gordon 
field equations (5) with the metric (9). In the Gowdy model, these functions are periodic in x. However, in the 
Schmidt model i£l and we need to impose asymptotic conditions for the fields in the limits x — > ±00. Here we will 
make use of the same conditions introduced in [7] 

T^TiW + O^- 1 ), N^N^ + Oix- 1 ), 

t -» T±(t) +0(x- 1 ), N x ^O{x~ 1 ), (10) 
4>-> 0(x- v ). 

Due to these fall-off conditions for the fields the boundary integral in the action (6) vanishes in the non-compact case 
(it is trivially zero in the compact one). 
Introducing the metric (9) in the action (6) it is straightforward to write it in the canonical form 

S[g, cj>] = J J ( jf (p 7 7 + Pt t + P^) dx - H[N, N x }^j dt, (11) 

where the dot denotes dt, p-y, p T , and p^, are the canonically conjugate momenta [with the fall-off conditions inherited 
from (10)], and H[N, N x ] is the Hamiltonian that can be written as the sum of the first class constraints C[N] and 
C X [N X ] 

H[N,N X ]=C[N]+C X [N X ] = [ NCdx+ [ N x C x dx. (12) 

Jx Jx 

Here 

C := ^-c--/ 2 (V - 7 V - (8G3) 2 Pr Pl + { ^pl + r0' 2 ) , 
C x := e"T/ 2 {prT 1 - 2 P ; +p. n ' + <f/ P4> ) , 

and we have denoted d x by a prime. In the following we will choose units such that 8G3 = 1 so that the previous 
constraints become 

C[N] = J iVc-^ 2 ^2r" - j't' - p T p 7 + ^-p 2 + T<f>' 2 ^ dx, 

C X [N X ] = f N x e-^ 2 (p T r' - 2 P ; + Pll ' + cf>' P<t> ) dx, 
Jx 

where N and N x are arbitrary functions. As the Hamiltonian is zero -regardless of X- we have no dynamics so in 
order to proceed further we will deparameterize the theory and introduce a suitable phase space variable to play the 
role of time. 



B. Deparameterization and Reduced Phase Space 

Deparameterization has been discussed in general in many places (see, for example, [25]); for the system under 
consideration we closely follow [6, 7] where the reader is directed for details. The fact that symmetry forces the area 
density of the group orbits r to have time-like gradient and the analogy of the present models with the Einstein-Rosen 
waves 2 suggests that we choose t as our time variable by imposing 

t' = and p 7 + 1 = 0. (13) 



2 For Einstein- Rosen waves the area density has a space-like gradient and the gauge fixing that simplifies the interpretation of the model 
leads to use it as a radial coordinate. 
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The first condition implies that r is constant on the spatial slices of the foliation whereas the second basically means 
that the same is true for f . Note that the last condition is essentially equivalent to p'^ — because in the non-compact 
case we have imposed a fall-off behavior for p 1 [7] forcing it to asymptotically approach — 1. These conditions imply 
that all the degrees of freedom but a single one -that will change from one slice to another and will become the time 
variable of our model- are eliminated. In the compact case one must use a gauge fixing of the type p y + p = for 
a constant p. This takes care of a global degree of freedom that must be considered (see [6]). As we are mostly 
interested in the field-theoretic aspects of quantization for these systems we will drop this here by using the condition 
p 7 + 1 = in both cases. At this point one can check that the gauge fixing conditions (13) are admissible by 
computing their Poisson brackets with the constraints [7]. One can also see that there is just one first class constraint, 
say C[N — e 7 / 2 ], which cannot be solved after fixing the gauge. This is found by solving for the lapse and shift that 
give a zero Poisson bracket of the Hamiltonian with the gauge fixing conditions. In fact the choice N = e 7 / 2 will be 
such that {t, C[N = e 7//2 ]} = 1 and, hence, we can identify r with the time parameter t of the system. 
In these models there is a quadratic conserved momentum 



/ p<t,(f>' dx 
Jx 



(14) 



that is forced to be zero in the compact case. As a consequence of this the reduced phase space for the Gowdy model 
is non-linear. 

In order to complete the characterization of the reduced phase space of the deparameterized theory we must solve 
the first class constraints together with the gauge-fixing conditions. When this is done we find that the action of the 
system can be completely expressed in terms of the field <f> and its canonically conjugate momentum p^ to give 
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Jty JX 



dxdt. (15) 



The field equations are equivalent to the Hamilton equations derived from the time-dependent Hamiltonian 

H(t) = ^(§+¥ 2 )dz. (16) 

Once these equations are solved -and the constraint has been taken into account for Gowdy- the full four dimensional 
metric can be built from their solutions by following in reverse the reduction process for the metric and using the fact 
that the 2+1 dimensional metric (9), defined in (0, oo) x X x S 1 , becomes 

gab = c 7 (-V a tV 6 t + V a xV b x) + t 2 V a aV b a 

where 7' = p<p(f>' . This metric displays a curvature singularity at t — 0. Finally it is important to remark [6, 7] that 
the reduced action (15) corresponds, precisely, to that of a massless scalar field (independent of the a coordinate) on 
a fixed background spacetime with a metric given by 

g ab = - V„t V b t + \7 a x \7 b x + t 2 \7 a a \7 b cr. (17) 

This metric is defined on a 3-manifold (0, 00) xXxS 1 and shows, again, the singular behavior of the model at t = 0. 



III. EVOLUTION OPERATORS 



This section is devoted to study the quantization of quadratic, time dependent, field Hamiltonians such as (16). 
The main difficulty to obtain evolution operators in closed form is due to their explicit time dependence that precludes 
us from writing them as exponentials of the Hamiltonian. A first point to notice is that some of the functions of time 
that appear in (16) -say the one in the p 2 term- can be eliminated by a simple redefinition of the time variable. Let 
us consider, then, the problem of quantizing a system with a time dependent Hamiltonian of the general form 

H(t) = J (^Tr 2 (x)+A(t)ir(x)p(x)+u} 2 (t)(p' 2 (x)^jdx. (18) 

Here ip(x) and ir(x) are a field and its canonically conjugate momentum [{f{x),n(y)} = S(x,y)} defined on X, that 
may be either the real line R (Schmidt) or the circle S 1 (Gowdy). The prime denotes the x-derivative. The time 
functions A(f) and u(t) will be determined by the particular model considered. We allow for cross terms -absent in 
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(16)- involving fields and momenta as they will appear in the discussion of the unitarity problem. We diagonalize the 
Hamiltonian (18) by writing 3 



<p(x) 



cos (kx — —^jip(k) dp(k), 



(19) 



where ip(k) and w(k) are also canonically conjugate [{ip(k),n(q)} = S(k,q) with a Dirac or Kronecker delta]. Here X 
denotes the real line for X = M and the integers Z when X = S 1 . The measure dp(k) simply refers to the fact that 
the previous integral either extends to R or becomes a sum over the integers respectively. In terms of them we have 



H(t) = \j. (ir 2 (k) + 2A(t)Tr(k)(p(k) + k 2 u; 2 (t)(p 2 (k)^ dp(k). 



(20) 



This somewhat unusual diagonalization has the advantage of decoupling the different modes right from the start 
and avoids the problems encountered in [6] with the non-diagonal form of the Hamiltonian in terms of creation and 
annihilation operators. As we can see this is a sum of time dependent uncoupled Hamiltonians that arc closely related 
to the harmonic oscillator with a time dependent frequency (even when the A(t)n(k)(p(k) term is present). In fact, 
denoting d t by a dot, the Hamilton equations derived from (20) lead to 



with 



(f(k,t) + Q 2 (k,t)ip(k,t) = 



Q.\k,t) :=k 2 uj 2 (t)- A 2 (t)- A(i), 



(21) 



that we will suppose to be positive in the following. 

To quantize the model we promote the field and momentum to formal algebraic objects ip(k) and Tt(k) satisfying 
the canonical commutation relations [<p(k), n(q)] = ihS(k, q) and symmetrize the cross term in (20). In the main body 
of the paper we take ft = 1 . This way we obtain 



H(t) 



1 



7r 2 (fc) + A(i) [it{k)(p(k) + <p{k)n(k)} + k 2 uj 2 (t)^ 2 (k)jdfi{k). 



(22) 



Explicit solutions for the evolution operator for problems similar to this have been known for a number of years 
(see, for example, [26, 27]). Using them (see appendix A) we can write down the formal quantum evolution operator 
U(t,to) as the product 



U(t,t ) =V(t,t Q )S(t,t Q )Tl(t,t ) 



with 



T>(t,t ) := cxp 



<S(Mo) : = exp 



2^ 



p(k,t ) p(k,t) 



log 



p(k,t ) p{k,t) 
p(k,t ) 



K{t,to) := exp I - - 



p(k,t) 



dr 



t P 2 (k,r) 



n{k) 



P 2 (k,t ) 



+ A(t)-A(t )\tp 2 {k)dp{k)), 
p(k,t Q ) 



(23) 



p{k,t ) 

7r(fc) 



A(t ) j m 

p{k,t ) 



p 2 {kM) 



n(k) 



P(k,t ) 

p(k,t ) 
p(k,t ) 



A(t ) m 



0{k) dp(k) , (24) 



A(t ) m 



dp(k) . (25) 



3 We purposely use the same letter to represent a field and its transform. 
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Here p(k,t) is any solution to the Ermakov-Pinney equation 4 

p(k, t) + Q 2 (k, t)p(k, t) = p- 3 (k, t) (26) 

with n 2 (k, t) defined in (21). It can be shown (see appendix B) that p(k,t) never vanishes and U(t, to) is independent 
of the choice of this solution. It must be said that the previous way to factorize U(t,to) is by no means unique 
although this form is specially adapted to the discussion of its unitarity in the next section. 

If we give a representation for the fields ip(k) and n(k) satisfying [ifi(k),rr(q)] = iS(k, q) and such that the exponents 
(multiplied by i) in (23,24,25) are self-adjoint, it is straightforward to check that Ufa, to) = I, (/(t^to)^ 1 = U{t,to)\ 
U(t2,to) = U(t2,ti)U(ti,to), and idtU(t,to) = H(t)U(t,to) where H(t) is the quantum Hamiltonian (22), i.e. U(t,to) 
is the quantum evolution operator of the system in this representation. The evolution of the field and momentum 
operators in the Heisenberg picture can then be computed in a straightforward way and shown to satisfy the classical 
equations of motion as expected (see appendix A) . 



IV. FOCK REPRESENTATIONS AND UNITARY IMPLEMENTABILITY OF TIME EVOLUTION 



A. General Considerations 



As a first application of the previously derived formula we will discuss here the problem of finding Fock space 
representations for the field and momenta leading to unitary evolution operators. This is important because of the 
known obstruction to the implcmcntability of time evolution for this model in some of the descriptions considered in 
the past. Our results in this regard strongly support the satisfactory resolution of this issue that appears in [19]. 

Let us suppose that we take a Fock space T and write (overbars denote complex conjugation) 

ip(k) = f(k)a k + f(k)a{, Tt(k)=g(k)a k +g(k)a{ (27) 

in terms of the annihilation and creation operators a k and a\ with [afc,a|] = 6(k,q). The functions / and g must 
satisfy 

f(k)g(k) - f(k)g{k) = i, (28) 

[so that [<p(k),n(q)] = iS(k,q)] but are otherwise arbitrary at this stage. We want to find / and g such that T>(t,to), 
S(t,to), 7Z(t,to) -with normal ordered exponents to prevent the appearance of infinite phases- are unitary and 
differentiable in t. An efficient procedure to discuss this issue has been put forward by Torre in [17] by using the 
theory of unitary implementation of canonical transformations. Let us consider a general quantum operator of the 
following type 

exp (i J [Ai(fc) : £ 2 (k) : +A 2 (fc) : ir 2 {k) : +A 3 (fc) : ip(k)n(k) + it(k)<p(k) :] dp(k)^j , (29) 

where Ai ; 2,3 are real functions of k. Notice that T>(t,to), S(t,to), and 7l(t,to) are particular cases of it with a 
parametric dependence on t and to- Introducing (27) in (29), the operator can be written as 

exp ^^[xi(fc)a fc a fc + xi(fc)aiai + 2 X2(fc)aiafc]d^(fc)^ , (30) 

where %i := ^i.f 2 + ^2.9 2 + 2A3/5, an d Xi '■— ^i|/| 2 + -Msl 2 + M(fS + fo) is a re& l function. We want to know 
now if the exponent (times i) defines a self-adjoint operator. For fixed values of t and to this is done by studying the 
auxiliary dynamics -in a fictitious time parameter s- defined by the exponent taken as a classical Hamiltonian 



/ lxi(k) 
Jx 



afeflfc + Xi(k)a k a k + 2x2(k)a k a k ] dp(k), (31) 



4 The Ermakov-Pinney equation is of central importance in the classical and quantum treatment of the harmonic oscillator with time- 
dcpcndcnt frequency. It has been used, for example, in some field theoretic cosmological problems (see for example [28, 29]) and 
minisuperspace models [30, 31]. In this paper we use it in the context of the Fock quantization of time dependent field theories. 
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and using the results on quantum implcmcntability that appear in [17]. The modes a k and a k in (31) are defined by 
the classical field and momentum ip(k) = f{k)a k + f(k)a k , n(k) = g(k)a k + g(k)a k and satisfy {a k ,a q } = —id(k,q). 
In practice it is convenient to consider the evolution equations 

^ = {a k ,F} = -2i (X2(k)a k + Xi(k)a k ) , ^ = {a k , F} = 2i {x2{k)a k + Xi(k)a k ) 

that are equivalent to the second order equations 

d 2 a k 



ds 2 



= 4(| X i(fc)| 2 - xi(k))a k = 4 (A 2 (fc) - X 1 (k)X 2 (k)) a k . (32) 



These are linear equations so their solutions have a linear dependence on the initial conditions a k (so) and da k /ds(so) = 
—2i (x2(k)a k (so) + xi(k)a k (so)) and, hence, on a k (so) and a k (so). In order to guarantee unitary implementability 
it suffices to show that the integral over X of the modulus squared of the coefficient of a k (so) that appears in the 
solution of (32) is convergent. Finally, in order to verify the strong continuity of the transformation in the auxiliary 
parameter s we have to check that the following limit 

lim / \a k (s) - a k {s )\ 2 dp(k) = (33) 

holds for the solution a k (s) of (32) with square summable initial data a k (so). We will obtain now general conditions 
that guarantee that the previously defined operators T>(t,t ), S(t 7 t ), and TZ(t,t n ) (23)-(25) are unitary: 

• T>(t,to) is a quantum operator of the form (29) with A 2 , 3 = and Ai ^ given in terms of the solution to the 
Ermakov-Pinney equation (26) by 

i ( p{k,t) \ i f p(Mo) A 

In this case the solution to (32) is given by 

a k (s) = [1 - 2i(s - s Q )x2(k)]a k (s ) - 2i(s - s a )xi(k)a k (s ) 
with xi = X\f 2 an d X2 = Ai|/| 2 . Then, the unitarity condition is 

/ Af(fc)|/(fc)| 4 d M (fc) < oo, (34) 

and the strong continuity condition (33) becomes 

= lim(.s- So ) 2 / A 2 (fc)|.f(fc)a fc ( S o) + |/(fc)| 2 a fc ( S o)| 2 d^(fc). (35) 

s^s a Jjr 

This is trivially satisfied whenever the last integral is well defined. 



S(t, to) is a quantum operator of the form (29) with A 2 = 0, Ai l3 ^ 0: 

Mk) = -(^-A(t )W (M0) 



\p(k,t ) ) p(k,t) ' 

\ n \ 1 i /°(Mo) 

As(fc) = 2 log 7(My 

The solution to (32) is now 

flfe(s) = {cosh[2A 3 (fc)(s - s )] - ix 2 (k)X^ 1 (k) sinh[2A 3 (fc)(s - s )}} a k (s ) 
- iXi(k)K sinh[2A 3 (fc)(s - s )] a k (s Q ) , 

where X i - XiP + 2X 3 fg, and X 2 = A x |./| 2 + X 3 (fg + fg). 
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The unitarity condition in this case is 

|/(fc)| 2 |2. 9 (fc) + Ai(fc)/(fc)Aa '(fc)! 2 sinh 2 [2A 3 (fc)( s - «„)] o>(fc) < oo. (36) 



It is important to point out here that the condition (28) implies that |/(fc)| 2 |g(fc)| 2 > 1/4 and, hence, the 
convergence of (36) requires an appropriate fall-off of |1 + Ai// (2A 3 g)| 2 sinh 2 [2A 3 (s — sq)]. The strong continuity 
condition is now 

= lim / |{cosh[2A 3 (fc)( S - s )] - i[f(k)g(k) + f(k)g(k) + A 1 (fc)A 3 - 1 (fc)|.f (fc)| 2 ] sinh[2A 3 (fc)( S - s )} - l}a k (s ) 
-*[A 1 (fc)A 3 - 1 (fc)/ 2 (fc) + 2f(k)g(k)} sinh[2A 3 (fc)( S - s )]a k (s )\ 2 d^(fc). (37) 

i Finally lt{t, to) is a quantum operator of the form (29) with Ai,2,3 ^ given by 

dr 



(k) P 2 (Mo) f* dr 

A 3 (fc) = f 4*1*4 - A(^ ^ 



t P 2 (fc,r) 



p(Mo) ' / 7t P 2 (fc,r) ' 



with A 2 (fc) := X 1 (k)X 2 (k)-Xt(k) = \^ t -^-^j >0. 



The solution to (32) is given by 

a k (s) = {cos[2A(fc)(s - s )] - ix2(fc)A -1 (fc)sin[2A(/c)(s - s )]}a fc (s ) 
- ixi(k)\^ 1 (k) sin[2A(fc)(s - s )] a fc (s ) , 

with xi = Ax/ 2 + A 2 . 9 2 + 2A 3 /<?, and X 2 = Axl.fl 2 + X 2 \g\ 2 + X 3 (fg + fg). 
Then, the unitarity condition is 

/ |A 1 (fc)/ 2 (fc) + A 2 (% 2 (fc) + 2A 3 (fc)/(%(fc)| 2 A- 2 (fc) sin 2 [2A(fc)(s - s )] dfi(k) < oo, (38) 
Ji 

and the strong continuity condition is 

= lim f |{[A 1 (fc)|/(fc)| 2 + A 2 (fc)|. 9 (fc)| 2 + A 3 (fc)(/(fc)g(fc) + /(%(fc))]sin[2A(fc)(s-s )] 

s^s„ j£ 

+iA(fc)(cos[2A(fc)(s - s )] - l)}a fe (s ) 
+ [A 1 (fc)/ 2 (fc) + \ 2 {k)f{k) + 2X 3 (k)f(k)g(k)} sin[2A(fc)(s - s )]a fc (s )| 2 \- 2 (k)dfi(k) . (39) 

B. The Gowdy and Schmidt Models 

As we have previously shown in section II, the reduced Hamiltonian (16) for the Gowdy and Schmidt models is 

H(t) = / x (^ + ^w)<i- 
After redefining the time parameter as t — e T and using the cosine transform (19) we get 

H(T) = \J (pl{k) + fc 2 e 2 V(fc)) d M (fc). (40) 
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In both cases there is a conserved quantity (14) of the form 

P = J kp4,(-k)<f>(k)dn(k). 

For the Gowdy model P is constrained to be zero. This is usually taken into account in the quantum theory as 
a condition on the physical states that is trivially preserved under the time evolution defined by the Hamiltonian 
(whenever it is well defined). 

We discuss now the unitarity of the time evolution defined by the explicit operator U(T, Tq) introduced above. The 
Hamiltonian (40) belongs to the general class (20) with A(T) = and uj 2 (T) = e 2T . Then (21) gives ft 2 (fc,T) = k 2 e 2T 
and, in this case, the general solution to the Ermakov-Pinney equation (26) can be written in terms of Bessel functions 
as 

p(k, T) = [AJ 2 (\k\e T ) + 2BJ (\k\e T )Y (\k\e T ) + CY*(\k\e T ) r ' 2 

with AC — B 2 = 7r 2 /4. It is possible to understand the difficulties to get a unitary evolution operator as the impossi- 
bility of fulfilling the unitarity conditions for some of the operators used to build U(T, T ), in particular the condition 
(36) for S(T,Tq). To show this we look at the leading asymptotic behavior of Xs(k) = (1/2) log[p(fc, 7b)/p(fc, T)] 
when \k\ — > oo that is given by (T — Xb)/4 plus a bounded periodic function of |fc|. We have also X\(k)/\3(k) — 
—2[p(k,To)/p(k,To) — A(To)] <~ 1. This shows that the convergence of (36) would require lim^^^ f(k)/g(k) = —2 

which is not compatible with (28) and, hence, S(T,T ) cannot be unitary. 

The way this factor fails to be unitary suggests that one could avoid this problem by introducing from the start 
(say at the Lagrangian level and, hence, before obtaining the Hamiltonian) a new field £ differing from in a certain 
time dependent factor. The reason is to modify the leading asymptotics of A3 (A;) to improve the convergence of (36). 
A possible way to change that behavior within the conceptual scheme that we are using here is to factorize the scalar 
field (f>(k, T) = h(T)£(k, T) for an appropriately chosen function h(T). It is always possible to find a unique (modulo a 
multiplicative constant) time redefinition for a given h in such a way that the new field £ satisfies, again, a differential 
equation for a harmonic oscillator with time (and k) dependent frequency. It could be possible, in principle, to obtain 
general conditions for h that guarantee that the function As(fc) has the right asymptotic behavior in k. Physically one 
would expect that a choice for which the frequency squared is a sum of k 2 plus a^ decreasing function of time would 
work as the system would approach a free one for which the evolution operator U in the form written above is well 
defined and unitary. The essentially unique way to do this (see appendix C) is to introduce a new field £ satisfying 

£(x, t) = Vi(j){x, t) (41) 

as was done in [14, 19]. Here the appropriate time variable is precisely the original one t — e T . The Hamiltonian in 
terms of £(x) and its canonically conjugate momentum p%(x) becomes [14] 

Notice that this is not the Hamiltonian considered in [19]. Although it is related to it by a canonical transformation, 
the two quantum dynamics are different because the same wave function at an initial time to evolves differently. 
Diagonalizing (42) by using the transform (19), defined now for £ andp^, and introducing the corresponding operators 
and Pz(k) we get 

H(t) = IJ_ (pi(k) + ^mrnik) + iwtw] + dp(k). 

As we can see the Hamiltonian (42) belongs to the class (18) with A(t) = l/(2i), to 2 {t) = 1. In this case, by using 
(21), n 2 (k,t) = k 2 + l/(4£ 2 ) and a solution to the corresponding Ermakov-Pinney equation (26) is 

p(k, t) = ^F/2[J 2 (\k\t) + Y 2 (\k\t)] 1 / 2 . (43) 



The asymptotic behavior of the functions associated to this solution that appear in the unitarity conditions (34,36,38) 
is given in the following table 
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Function 



1*1 



oo 



\k\^0 



p(k,t) 



it 

-log(|fc|t) 

IT 



dr 



t P 2 (k,r) 



C(t,t )\k\ 



log(*Ao) 
log 2 |fc| 



log 



p(Mo) 
p{k,t) 



16\k\W 



log(toA) 



1 1 

2 + \^i\k\_ 



p(M) 
P {k,t) 



-A(t) 



1 



1 



8|/c| 2 i 3 2* 



*Iog(|fc|t) 



where we do not need the specific form of C(t,to). They can be easily derived from (43) with the exception of the 
second one that requires the use of equation (B2) of appendix B. Choosing, for example, 

/(*;) = l/V2(l + |fc|) and g(k) = -iy/(l + \k\)/2 

it is straightforward to check that the unitarity conditions (34,36,38) are satisfied together with (28). Finally, strong 
continuity in the auxiliary parameter s follows from (35,37,39), by using the asymptotic expansions given above and 
the fact that afc(so) is square summable. We conclude, hence, that the exponents (times i) in T>(t,to), S(t,to), and 
TZ(t,t a ) are self-adjoint and these are unitary 

In order to verify that U(t, to) satisfies the evolution equation one must be able to compute derivatives in t. This 
can be shown to be formally possible due to the fact that t derivatives of the exponent trivially commute with it. 
Also, one can see (appendix A) that the evolution of the field and momentum operators generated by U(t,to) is 
diffcrentiablc and their time derivatives are obtained as commutators with the Hamiltonian. This allows us to make 
sense of the time derivative of U(t,to) indirectly through its action on the basic operators. Of course, one can try to 
study the behavior in (t, to) in a rigorous mathematical sense from first principles but this is beyond the scope of the 
present paper. 



V. CONCLUSIONS AND COMMENTS 



Some comments are in order now. First of all it must be said that if one is only interested in Fock space representa- 
tions the evolution of n-particle states can be simply derived from the evolution of creation operators in the Heisenberg 
picture (that can be directly read from the classical dynamics of the system) and the evolution of the vacuum state. 
Furthermore, the vacuum evolution can be written in closed form as in [32, 33]. This provides the matrix elements 
of the evolution operator in a concrete Hilbert space basis. A potential advantage of the framework discussed in this 
paper is that one can, in principle, try to use non-Fock representations to define the evolution operator given in closed 
form by (23,24,25). A second comment is that 

^2 = -±= j f [f(k)a k + f(k)a{] cos (kx J) d M (fc), 

interpreted as an explicitly time dependent operator in the Schrodinger picture evolves, in the Heisenberg picture, 
as the classical scalar field <fi(x, t) that encodes the physical degrees of freedom of the Gowdy and Schmidt models. 
This evolution is perfectly well defined and unitary. As we can see the problem is not that it is impossible to define 
a unitary evolution for an object that behaves as the scalar field <j> but rather the impossibility of doing this with a 
representation of the type 

0(fc) = f(k)a k + f{k)a\ , p^{k) = g(k)a k + g{k)a\ . 

We want to emphasize that the use of the field £ defined in (41) is as justified as that of <j> because it can be chosen 
as the fundamental object already at the Lagrangian level. It is important, however, to realize that the quantum 
dynamics must then be obtained from the corresponding Hamiltonian. In this respect it should be noticed that it is 
possible to write down the Hamiltonian corresponding to the dynamics considered in [19], build the unitary evolution 
operator in that case with our methods. 

The possible uses of the approach that we have presented here are manifold. In addition to the discussion of 
consistency issues related to the unitarity of the time evolution of the system other problems that may be tackled 
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are the generalization of the coherent and squeezed states for the harmonic oscillator with time dependent frequency 
to these quantum cosmological models, the study of particle creation, and the introduction and discussion of new 
representations to quantize the system. 
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APPENDIX A: COMPUTING THE UNITARY EVOLUTION OPERATOR 



We give here some details about the construction of the formal evolution operators for the kind of system considered 
in the paper. In order to quantize the model defined by the Hamiltonian (20), we promote the field and momentum to 
formal operators <p(k), 7r(fc), satisfying the canonical commutation relations [ip(k),jr(q)] = iH5(k,q), and symmetrize 
the cross term in the classical Hamiltonian to arrive at (22). We want to solve now the Schrodinger equation 5 
ihd t \ip{t)} — H{t)\^p{t)). The strategy that we follow is to generalize the results already known for a single harmonic 
oscillator with a time dependent frequency to an infinite system of uncoupled harmonic oscillators [26, 27]. In order 
to obtain the quantum evolution operator, we define 

D(t; f) := exp (-± jf f(k, t)£ 2 (k) d/x(fc) 

S{t;g) := exp (j^ ^g{k,t)[£(k)^k) + n(k)<p(k)} d,i(fc) 

depending on some functions / and g that will be fixed in the following. Let us consider their product T(t; f, g) = 
S(t; g)D(t; /) and introduce the state vector \x(t',f,g)} — T(t; f, g)\ip(t)) written in terms of a solution to the 
Schrodinger equation \ip(t)}. We find now the equation that |x(i; /,£/)) must satisfy and choose the functions f(k,t) 
and g(k,t) in such a way that the evolution for this vector can be easily obtained 



iWt \ X (t; /, g)) = [f{t; f, g)H{t)T* (t; /, g) - ihf(t; /, g)d t f1 (t; /, g)j \ x (t; /, g)) 

= dn(k)(c- 29( - k ^Tr 2 (k) + f(k,t)+A(t)-g(k,t) <p(k)it(k) + n(k)0(k) 
'.f(k, t) + fcV(t) + 2A(t)/(fc, t) + f(k, t)} ^(k)) | x (i; /, .9)) 



2g(k,t) 



+e 



where the dot denotes d t and we have used the following relations satisfied by D(t; /) and S(t;g) 

D(t;f)<p(k)Dl(t;f) = <p(k), 
D(t; /)7r(fe)Dt(t; /) = TT(fc) + f(k, t)0(k), 
S{t;g)<p{k)&(t;g) = <p{k)e 9 ( k ' t \ 
S(t;g)n(k)SHt;9) = me- 9ik ' t) - 
The form of (Al) suggests to choose functions /, g such that 

f(k,t)+A(t)=g(k,t), 

/ 2 (fc,t) + kW(t) + 2A(t)/(M) + /(*,*) = e- 4 ^. 



(Al) 



5 We assume that a suitable Hilbert space representation exists. 
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The equation for \x(t)} simplifies then to 

ihd t \ X (t)) = \l d/x(*) e~ 2 ^ [% 2 (k) + <p 2 (k)] \ X (t)). 

A concrete way to achieve this is to take g(k,t) — logp(fc,i) [and, hence, f(k,t) = p(k,t)/ p(k,t) — A(t)] with the new 
function p(k, t) obeying the auxiliary Ermakov-Pinney equation 

p(k, t) + Q 2 (k, t)p(k, t) = p- 3 (k, t). (A2) 

Here the dot denotes dt and Vi 2 (k,t) = k 2 u> 2 (t) — A 2 (t) — A(t), that we will suppose to be positive in the following. 
As we can see this equation is closely related to the classical equation of motion for a harmonic oscillator with time- 
dependent frequency f2(fc,i). In terms of this new function the state vector \x(t)) satisfies the differential equation 



mm) = \ ( K J ± -^j ) [^ 2 (k)+'P 2 (k)] d M (*o) \ X (t)) 

that can be easily solved to get 

\ X (t)) := to)|x(to)) = cxp JL Qf* ^7)) [* 2 « + ^W] d ^)) lx(*o)>- 

At this point one can build the solution to the original Schrodinger equation from \x(t)) an d the operator T(t). 

The solutions to (A2) have some interesting features, in particular they never vanish and are such that |x(t)) is 
independent of the initial conditions chosen for p(k, t), as will be shown in appendix B. Other choices for the function 
g(k,t) allow the formal solution of the equation for |x(0) but are problematic because the evolution operators thus 
obtained are not well defined for all values of t [27]. Finally, going back to the original state vector \tp(t)), we obtain 
the formal quantum evolution operator 6 

U(t,t ) =f\t)R(t,to)f(t ), 

with f(t) = S(t)D(t) and 

%)=cxp (ri/i 

S(t) = exp QL j( logp(M) [<p(k)n(k)+n(k)<p(k)} d^k^j , 

R(t,t o) = exp (-±1 (jf -Jjl_) [**(*) + M k)) . 

It is clear, by construction, that if we can find a suitable representation for <p(k) and Tt(k) as self-adjoint operators on 
a certain Hilbert space, U(t,to) satisfies 

U(t ,t ) =1, 

{/(Mo)" 1 = U(tM)\ 
U(t 2 ,t ) = U{t 2 M)U(tut ), 
iW t U(t,t ) = H{t)U(t,t ) 

where H(t) is the quantum Hamiltonian (22). It is interesting to notice (see appendix B) that 

\Jt P 2 (k,r)J p{k,t)p{k,t ) 



p(k,t) 



A(t) 



ip 2 (k) d/i(fc) 



The evolution of a state vector from time to to t is given by \ip(t)) = U(t, to)IV'(to))- 
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where u(k,t) is the unique solution to the homogeneous equation u[k,t) + Vl 2 {k,t)u{k,t) = satisfying u(k,to) = 
and ii(k,to) = 1. This allows us to write explicit and closed expressions for U(t,t n ) whenever we can give closed 
solutions to the previous equation for u(k, t) because solutions to the Ermakov equation can be obtained from these. 

In order to discuss the conditions that the representations for the field and momentum must satisfy it is convenient 
to factorize it as U(t,t ) = T>{t,to)S(t,to)'R,{t,tQ) with 

V(t,to) :=Dt(t)D(t ), 

n(t,t Q ) := D\to)&(to)R(t,to)S{to)D(to). 
These expressions directly lead to (23,24,25). The evolution of the field is 

<p H {k,t ,t) = p(k,t) 



1 ( f* dr \ ,, ./p(k,t ) A . A . / I 1 * dr 

" P(Mo) Trn-A(fo) sin' ' 



_p(k,t ) \Jt P 2 {k,T)J ' \p(k,t ) J \J to p 2 (fc,r) 

dr 

Similarly we can compute the momentum operator in the Hcisenbcrg picture 



+ p(k,t)p(k,t Q ) sin [ I _^ ) n(k). (A3) 



^-(k,t ,t) = ~[<p H {k,to,t),H H (t)] =n H (k,t ,t). 

It is straightforward to show that these are formal solutions to the classic equations of motion -as one would expect 
given the quadratic nature of the Hamiltonian- because the coefficients of <p(k) and n(k) are solutions to the classical 
equations <p(k,t) + Q 2 (k,t)(p(k,t) = 0. 



APPENDIX B: SOME PROPERTIES OF THE ERMAKOV-PINNEY EQUATION 

We discuss in this appendix some properties of the Ermakov-Pinney equation 

p(t) + n 2 (t) P (t) = p- 3 (t) 

that are relevant to show the independence of the evolution operators from the choice of initial conditions. The most 
general analytic solution to this equation 7 can be written as 

p{t) = + [Au 2 (t) + 2Bu(t)v(t) + Cv 2 {t)\ 1/2 

where u(t) and v(t) are linearly independent solutions to the time-dependent harmonic oscillator equation ip{t) + 
^l 2 {t)ip{t) = 0, A,B,C are constants satisfying AC — B 2 = 1/W 2 [u,v], and W[u, v] :— uv — uv denotes the (time 
independent) Wronskian of the pair of solutions u and v. The square of a solution can be expressed as 



p 2 (t) = [«(t) v(t)] 



A B 




u(t) 


B C 




v(t) 



where (u, v) ^ (0, 0) for all t (zeros of solutions u(t) and v(t) alternate and never coincide as a consequence of Sturm's 
theorem). Hence 

= AC - B 2 = W~ 2 > 0, 
therefore, the quadratic form p 2 {t) is positive definite and p(t) > for all t. 



A B 
B C 



Here we drop the parameter k that appears in (A2) because the results are trivially generalized to that case. 
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The general solution of the time-dependent harmonic oscillator equation tp(t) + Q 2 (t)<p(t) = can be expressed as 



<p(t) 



Po 



cos 



to 



dr 
7F) 



p(t)p sin ^ 



dr 
~pHj) 



ipa + p(t)p sin 



to 



dr 



(Bl) 



where pit) is an arbitrary solution to the associated Ermakov-Pinney equation p(t) + tt 2 (t)p(t) = p~ 3 (t); the subindcx 
denotes evaluation at to. Notice that the positivity of p{t) guarantees that the previous expression is well-defined 
for any value of t. As ip(t) is uniquely determined by the initial conditions ipo and ipo the functions 



cos 



to 



dr 
~pHr) 



must be independent of the choice of the solution p(t) to the Ermakov equation. This result has some interesting 
consequences regarding the quantum evolution operators studied in the paper. First of all it is possible to show that 
the unitary evolution operator for a single quantum harmonic oscillator with time-dependent frequency (obtained 
with the techniques explained in the paper) is independent of the choice of p(t) . This is obvious if one looks at the 
evolution for the field operator given by (A3). Second, it is easy to see from (Bl) that the solution u(t) to the equation 
u{t) + Q 2 (t)u(t) = satisfying u = and iio = 1 is 



u(t) = p(t)p sin 



dr 

toTF) 



and, hence, 



sin 



to 



dr \ u(t) 

7^))~pWo' 



(B2) 



As the integrand in the previous expression is positive, the integral is a strictly increasing function of t. The value of 
the integral itself is given in this case by the unique continuous and increasing extension to M of the arcsin function 
with argument u(t)/[p(t)p ]. This is used to obtain one of the asymptotic expansions appearing in IV B. 



APPENDIX C: ON THE DEFINITION OF THE FIELD £ 

We will discuss in this appendix the obtention of the field £ introduced in section IV B. The reduced phase space 
for the Gowdy and Schmidt models can be coordinatized by (4>(k) , p<p(k)) and the Hamilton equations derived from 
(40) lead to the time-dependent harmonic oscillator equation 

d%cj)(k,T) + k 2 e 2T <P(k,T) = 0. (CI) 

It is always possible to write (CI) in terms of a new field £ defined by 

4>(k,T) = h(T)£(k,T) , with h>0, 

to give 

dU(k, T) + 2^d T ak,T)+(^l + k 2 e 2T ^ * (*, r) = ; (C2) 

here h! and h" denote the first and second derivative of the function h respectively. After performing a time redefinition 
i(T), with t'(T) > 0, the equation (C2) can be written as 

d 2 i{k, t) + t'- 2 (T) (t"(T) + 2^P) t) + t'- 2 {T) + k 2 e 2 ^ £(fc, t) = . (C3) 

Although (C3) is not yet a time-dependent harmonic oscillator equation, it is always possible to eliminate the (<9 t £)- 
term by imposing 

t'h' 

t" + 2-^=0. (C4) 
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The general solution of (C4) is 



t(T)=to£j^- y (C5) 



'T 

where to > and T are constants of integration. With these assumptions the equation (C3) becomes 

<fr(M) 1 ( k " {T) 



{t>i(T)h(T) +k wfj)^- ' 

with t(T) given by (C5). This last equation is now of the form d 2 £(k,t) + £l 2 (k,t)£(k,t) = 0, with 

r> 2 (M)= JCiIL + fc2 ' " 



t' 2 {T)h{T) t' 2 {T) ' 

for any choice of h. Finally, it is possible to fix the function h by requiring that f2 2 is a sum of k 2 plus a function of 
time. Explicitly, by demanding e 2T /t' 2 (T) = 1, we get 

t(T) = e T ~ e T ° , h(T) = Vk>e- T / 2 (C6) 

or, equivalcntly, 

Particularizing (C6) for i = 1 an d Tq — > — oo we finally arrive at (41) 

f(M) = Vi<t>{k,t), 

where the field £ satisfies 

9^(fc,0+(^2+fc 2 V(fc,*)=0. 
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